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Abstract 

Recently, it was observed that self-interacting scalar quantum field theories 
having a non-Hermitian interaction term of the form g(i(f>) 2+s , where 5 is a 
real positive parameter, are physically acceptable in the sense that the energy 
spectrum is real and bounded below. Such theories possess PT invariance, 
but they are not symmetric under parity reflection or time reversal separately. 
This broken parity symmetry is manifested in a nonzero value for (</>), even if 
5 is an even integer. This paper extends this idea to a two-dimensional super- 
symmetric quantum field theory whose superpotential is S(4>) = —ig(icf)) 1+s . 
The resulting quantum field theory exhibits a broken parity symmetry for all 
5 > 0. However, supersymmetry remains unbroken, which is verified by show- 
ing that the ground-state energy density vanishes and that the fermion-boson 
mass ratio is unity. 
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I. INTRODUCTION 



It is extremely difficult to perform conventional perturbative calculations in supersym- 
metric field theory because expansions in powers of the coupling constant are infrared diver- 
gent. Furthermore, introducing a regulator in the form of a momentum cutoff or a lattice 
spacing to control this divergence breaks the supersymmetry invariance. An especially sim- 
ple way to solve this problem is to use the delta expansion, a perturbative expansion in 
powers of the degree of the nonlinearity of the interaction term [|IJ. 

The key idea of the delta expansion is to replace a self-interaction term of the form gcf) 4 
by g((f) 2 ) 1+6 , where 5 is regarded as a small parameter. Thus, the parameter 5 measures 
the departure from linearity of the field equation. A graphical procedure for expanding 
the Green's functions of a quantum field theory in powers of 5 is given in Ref. 0. The 
advantage of such an expansion is that it has a nonzero radius of convergence and that it 
yields numerically accurate nonperturbative information when 5 is set equal to 1. 

The delta expansion is broadly useful for nonlinear problems. It has been applied to 
many classical nonlinear ordinary differential equations of mathematical physics and it has 
given superbly accurate numerical results 0]. The delta expansion has also been successfully 
applied to some nonlinear partial differential equations of mathematical physics ||. 

In the context of quantum field theory the delta expansion has been used to study 
renormalization ||], local gauge invariance ||, stochastic quantization ||, finite-temperature 
field theory 0, and the Ising limit of quantum field theory 0. 

The delta expansion is particularly well suited for studying supersymmetric models of 
quantum field theory because expansions in powers of 5 are not infrared divergent and the 
delta expansion respects the supersymmetry exactly. Calculations to second order in powers 
of S have been done for the ground-state energy || and the fermion-boson mass ratio [T(| . 



Recently, we have examined a new class of scalar quantum field theories that are not 



symmetric under parity reflection |TT|| . The Euclidean-space Lagrangian for this class of 
theories is 

£ = ^(<90) 2 + \m 2 <? - g^f^ (5 > -2). (1.1) 

The Hamiltonian for such theories is not Hermitian. However, there is strong evidence that 
these theories possess energy spectra that are real and bounded below.[] The theories in 



There is copious numerical evidence that such theories possess positive spectra when m ^ 0. 
Furthermore, one can understand positivity from a theoretical point of view. Consider, for example, 
the weak-coupling expansion for the d-dimensional Euclidean quantum field theory defined by the 
Lagrangian C = ^(dcf)) 2 + ^m 2 (j) 2 +gi(f> 3 . For a conventional gcf) 3 theory the weak-coupling expansion 
is real, and (apart from a possible overall factor of g) the Green's functions are formal power series 
in g 2 . These series are not Borel summable because they do not alternate in sign. Nonsummability 
reflects the fact that the spectrum of the underlying theory is not bounded below. However, 
when we replace g by ig, the perturbation series remains real but now alternates in sign. Thus, 
the perturbation series is now summable and this suggests that the underlying theory has a real 
positive spectrum. 
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Eq. ( |1.1|) are a natural field theoretic generalization of a remarkable quantum mechanical 



Hamiltonian studied by D. Bessis and C. Itzykson: |12] 

H=^p 2 + ix 3 . (1.2) 

The Lagrangian in Eq. (|LlD is intri guing because it is not parity symmetric. This is man- 
ifested by a nonzero value of {(f)). It is interesting that this broken symmetry persists even 



when 5 is an even integer [fLl|j 



That fact that the Lagrangian in Eq. (1.1) has a broken parity symmetry suggests that 



if we construct a two-dimensional supersymmetric quantum field theory by using a super- 
potential of the form 



5(0) = -ig(i<P) 1+ \ (1.3) 

the resulting theory will also have a broken parity symmetry. The supersymmetric La- 
grangian resulting from the superpotential (|1.3|) is 



c = l -{dct>f + + l -S'{4>)^ + \[S{ct>)] 2 
= ^0) 2 + \$W> + \g{i + S)W)'W> - \g 2 (i<t>) 2+2S , (i-4) 

where ip is a Majorana spinor. 

The Lagrangian (|1.4j) raises an interesting question. Will the breaking of parity symmetry 
induce a breaking of supersymmetry? To answer this question we calculate both the ground- 
state energy -E gr0 und state and the fermion-boson mass ratio R as series in powers of the 
parameter 5. We find that through second order in 5, -E grcm nd state = and R = 1, which 
strongly suggests that supersymmetry remains unbroken. Based on the experience with 
these calculations we believe that our results are valid to all orders in powers of 5. It is quite 
difficult to break supersymmetry |13| 



This paper is organized very simply. In Sec. |l|we explain our calculational procedure. We 
derive a set of Feynman rules for obtaining the Green's functions to a given order in powers 
of 5 for the scalar field theory Lagrangian Next, in Sec. |T|we apply the procedures of 



Sec. to calculate the ground-state energy and the fermion-boson mass ratio to first order 



in 5 for the supersymmetric Lagrangian ( |1.4|) . In Sec. [IV] we perform these calculations to 
second order in 5. Finally, in Sec. [V] we examine some of the formal cancellations that occur 
in our calculations to see if there could be anomalous contributions to the ground-state 
energy. 



II. DELTA EXPANSION FOR A PARITY- VIOLATING SCALAR FIELD THEORY 

In this section we explain how to calculate the Green's functions for the scalar quantum 
field theory described by the Lagrangian C in Eq. ( |1 . 1| ) . Specifically, we follow the proce- 
dures described in Ref. [0 and determine a set of graphical rules for constructing the delta 
expansion through second order in powers of 5. These rules consist of the amplitudes for 
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the vertices and lines of a related Lagrangian in which the boson fields are raised to integer 
powers. 

We begin by expanding the Lagrangian £ in Eq. ( |1 . 1| ) to second order in powers of 5: 



£ = \(d<P) 2 + \m 2 <f + g<f + 5g<p 2 ln(z0) + ~5 2 # 2 [ln(i0)] 2 + 0(5 3 ). 
Observe that at 5 = the Lagrangian contains the coupling g: 

C Q = \(d<P) 2 + ^m 2 <p 2 + g<p 2 . 



(2.1) 



(2.2) 



Thus, an expansion in powers of 5 is clearly nonperturbative in g. 

In general, the functional integral representation for an n-point correlation function has 
the form (modulo the usual normalization factor) 



(0\<l>(x 1 )<f>(x 2 )<l>(x 3 )...<l>(x n )\0) = { D4 > 4 > (x 1 )4 > (x 2 )4 > (x 3 )...4 > (x n )e-I dxC . 



We expand the exponential in the integrand in Eq. (|2.3| ) and obtain 

e -JdxC = e -JdxCo S ^^ 

where 

£[(j)\ — 1 — 8g J dx (p 2 ln(i0) — -^> 2 9 J dx 2 [ln(i0 



(2.3) 



(2.4) 



dx<p 2 ln(i<; 



0(S 3 ). 



Next, we use the identity 



(2.5) 



(2.6) 



lnO'0) = ln(|0|) + -m- 



where (j)j\4>\ represents the algebraic sign of 0, to decompose the expression for 8 in ( [2.5| 
into its real and imaginary parts: 



£\ 



1 — Sg J dxcj) 2 ln(|0|) — -m5g J dx 4>\4>\ — -5 2 g J <ix0 2 [ln 
- h^5 2 g J dx<t>\(j)\ ln(|0|) + -7r 2 S 2 g J dx<f) 2 + \b 2 g 
+ ~iir5 2 g 2 f rfx0 2 ln(|0|) J dx <P\<P\ - ^n 2 5 2 g 



dx (J) 2 In 

2 



0(5-). 



(2.7) 



It is easy to construct an effective Lagrangian L having polynomial interaction terms 



that can be used to calculate the Green's functions of ( [2.7|) to first order in 5: 



C = C + \bg^+ 2 + Uv5ga<p 2 ^ + \ 



(2.8) 



To use this Lagrangian we must assume that the parameters a and 7 are integer. We can 
then read off a set of Feynman amplitudes: 



boson line: 



1 



p2 _|_ m 2 _|_ 2g ■ 



(2a + 2) boson vertex: --5g(2a + 2)!, 
(2a + 27 + I) boson vertex: — -i7i5ga(2a + 27 + 1)! 



(2.9) 



Note that the vertices are of order 5. Thus, if we are calculating to first order in S we need 
only include graphs having one vertex. Once a graphical calculation has been completed we 
then apply the derivative operator 



V 



d_ 
da' 



(2.10) 



followed by setting a = and 7 = 5- The technique of using a derivative operator to recover 
the Green's functions for the Lagrangian ( |1.1| ) is like the replica trick used in calculations 
in statistical mechanical models. It is a standard procedure used in all papers on the delta 
expansion and is is discussed in great detail in Refs. [^p|,|T0f . 

To perform calculations to second order in 5 we seek a higher-order effective Lagrangian 
£ of the form 



C = Co + (SA 1 + 5 2 A 2 )g ( j ) 2a+2 + (SB 1 + 5 2 B 2 )g<fi 2f3+2 

+ {Sd + b 2 C 2 )g^ 2 ^ + (5D 1 + 5 2 D 2 )g^ +2 ^\ 



(2.11) 



To determine the coefficients A\, A 2 , Bi, B 2 , C\, C 2 , D\, and D 2 , we replace £ by £ on the 
left side of Eq. ( |2.4| ) and expand to second order in 5. We then apply the derivative operator 

BI 



2 \ da dpi A \da 2 dp 2 



(2.12) 



and set a = 0, (3 — 0, and 7 = ~. By comparing with the right side of Eq. ( |2.7|) we obtain 
a set of seven simultaneous equations for the coefficients: 



T>(AtJ dx^ + B.J dx^ +2 
V (a 2 J dx <J) 2a+2 + B 2 Jdx <j) 2/3+2 

V (c, J dx 2 -+ 2 ^+ 1 + D 1 Jdx <p 2 ^+ l 

V (c 2 J dx 2a +^ +1 + D 2 J dx <j>W+*i+i 



<* = /3=0 
1=h 



1=h 



t* = /3=0 



rfx0 2 ln(|^)|), 



— ill I diXlf,y.'\ ■ 



-Z7T / dx(j)\(j)\\n(\(f)\), 
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V^A 1 Jdx<j) 2a+2 + B 1 Jdx<i) 2 P+ 2 j ^ = |cfe0 2 ln(|0|) 
vfc.J dx^^ + D.J dx^ +2 ^ ^ o = -ivr 2 ^ J dx^J, 

V^dJ dx<j) 2a + 2 J dx<f ) 2a+2 ^ 1 + A 1 D 1 J dx<p 2a+2 J dx<p 2 ^ 2 
+B l C l J dx^ +2 J dx<f>^ +1 + B 1 D 1 f dx^ +2 J dx<? 

= ^iir J dx (p 2 \n(\(j)\) J dx<p\<p\. (2.13) 



h2 7 +l 



^2/3+27+1 j 



*=/3=0 



The solution to these equations is: 



1 

2' 



1 

'2' 



Ci = —ilia, 
2 

J ° 1 = ~2 i7T ^ 



A 2 
B 2 



1 1 2 2 

4-8™' 



4 8 
r , 1. 

C 2 = -ITT, 

D 2 = -^tt- 



^2,o2 
7T /* , 



(2.14) 



We thus read off the Feynman amplitudes for the effective Lagrangian £ to second order 
in 5: 

1 

boson line: 



p2 _|_ m 2 _|_ 2g ' 



(2n -- 2) in. sun vertex: | --5 - + iyrVa 2 ) r/( 2n : 2)!. 



(2/5 + 2) boson vertex: (-5 - -5 2 + -n 2 5 2 (3 2 ) g(2(3 + 2)!, 

\2 4 8 / 

(2a + 27 + I) boson vertex: 



1 



1 



-in5a - -mb 2 1 #(2a + 27 + 1)!, 
2 4 / 

(2(3 + 27 + I) boson vertex: ( -m5j3 + -md 2 ] g(2(3 + 27 + 1)!. 



(2.15) 
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To second order in S one must include all graphs containing up to two vertices and treat the 
parameters a, (3, and 7 as integers. At the end of the calculation one must then apply the 
derivative operator V in Eq. ( j2.12| ) and seta = 0,/? = 0,7=|. 

In the next two sections we generalize this approach to the case of supersymmetric 
Lagrangians, identify the Feynman rules for calculating Green's functions, and use these 
rules to calculate the ground-state energy, the one-point Green's function, and the fermion- 
boson mass ratio. 



III. FIRST-ORDER CALCULATIONS FOR THE SUPERSYMMETRIC 

LAGRANGIAN 

In this section we show how to calculate to first order in 5 the ground-state energy 
and the fermion-boson mass ratio for the two-dimensional Euclidean field theory defined by 
Eq. (|1.4j). We follow the procedure described in Sec. [TJ and obtain a set of graphical rules 
for constructing the delta expansion to a given order in powers of 5. These rules are the 
amplitudes for vertices and lines of a related effective Lagrangian in which the boson fields 
are raised to integer powers. 

Our first objective is to obtain this related Lagrangian. We begin by expanding Eq. ( |1.4| ) 
to first order in 5: 

£ = \(d<P) 2 + l -g 2 ^ + ~z#V + \gH> 

+5g 2 <p 2 \n{i(j)) + -Sgijnp + ^Sg^ip ln(z</>) + 0(5 2 ). (3.1) 

We can display the real and imaginary parts of this Lagrangian explicitly by using the 
identity Q. 

Unfortunately, the Lagrangian (|3.1|) is nonpolynomial so we cannot read off a conven- 
tional set of Feynman rules. However, consider the following effective Lagrangian whose 
interaction terms are polynomial in form: 

C = ^(90) 2 + l -g 2 <P 2 + ~z#^ + l -g^ 

+ -5g 2 (j) 2+2a + -5(2a + l)g^(p 2a + -in6g 2 a0 2a+w + ^mbga^^ 2 ^ , (3.2) 
2 4 2 4 

where a and 7 are to be regarded temporarily as positive integers. Note that we recover the 
Lagrangian C in Eq. Q3.1|) if we apply the derivative operator V in Eq. fl2.10|) to C and set 
a = and 7 = \. 

To calculate the Green's functions for £ we use the same device. That is, we find the 
Green's functions for C. We then apply V to these Green's functions and set a = and 
7=2' ^° ca l cu l a te the Green's functions for £ we read off the following Feynman rules 
from Eq. Q: 

boson line: A(p) 



V 2 + 9" 
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(2a + 2) boson vertex 



fermion line 




--5g 2 (2a + 2)\, 
~5g(2a + l)\, 



(2a) boson and 2 fermion vertex 



(2a + 27 — 1) boson and 2 fermion vertex 



(2a + 27 + I) boson vertex 



--iix5g 2 a(2a + 27 + 1)! 
— m5ga(2a + 27 — 1)!. 



(3.3) 



Note that these rules are nonperturbative in the coupling constant g\ the parameter g appears 
nontrivially in all the vertex and line amplitudes. 

In addition to these Feynman rules one must remember to associate a factor of —1 with 
every fermion loop, and that in two-dimensional space gamma matrices are two dimensional 
so that the trace of the unit matrix introduces a factor of two. Moreover, one must be careful 
to associate with each graph the appropriate symmetry number. The fermions for this 
quantum field theory are Majorana fermions, which are nondirectional. Thus, for example, 
the symmetry number of a fermion self-loop consisting of a single fermion propagator is |. 

Using these graphical rules, we calculate below the zero-point, one-point, and two-point 
Green's functions to order 5. Because we are calculating to first order in 5 and all vertices 
in Eq. ( |3.3j ) are of order 5, we need only consider single-vertex graphs. We will show that 
to order 5 the ground-state energy vanishes and that the renormalized fermion and boson 
masses are equal. This indicates that the theory is supersymmetric to this order. We also 
see that the one-point Green's function is nonvanishing, which shows that parity is broken 
to order 5. These results are obtained using a formal manipulation of divergent quantities 
in the form of logarithmically divergent integrals. In Sec. [V] we re-examine these formal 
manipulations more carefully to see whether an anomalous structure could be present. 



The ground-state energy of a quantum field theory is the negative sum of the connected 
graphs having no external lines. Because we are considering graphs with no external lines, 
only the vertices in Eq. ( |3.3|) having an even number of boson lines contribute to the ground- 
state energy. Our results can be expressed in terms of a single dimensionless divergent 
integral A, which represents the amplitude of a boson self-loop (a loop formed from one 
boson propagator): 



where A(x) is the boson propagator in coordinate space [the Fourier transform of A(p)]. 
The amplitude for a fermion self-loop can also be expressed in terms of A: 



A. Calculation of the ground-state energy 




(3.4) 



(2*) 



1 




(3.5) 
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Two graphs of order 5 contribute to the ground-state energy density. These graphs are 
shown in Fig. |l[ The amplitudes for these graphs are the product of the vertex and line 
amplitudes multiplied by the appropriate symmetry numbers. The first graph is a pure boson 
graph consisting of one vertex with a + 1 boson self-loops attached. The amplitude for the 
vertex is — ^5g 2 (2a + 2)\, the symmetry number for this graph is 2~ Q ~ 1 / (a + l)!, the Feynman 
integral for the graph is A a+1 , and we include a factor of —1 because we are calculating the 
ground-state energy density. The second graph is a mixed boson-fermion graph consisting 
of one vertex with a boson self-loops and one fermion self-loop attached. The amplitude for 
the vertex is — ^5g(2a + 1)!, the symmetry number for the graph is 2~ a ~ 1 /a\, the Feynman 
integral for the graph is 2gA a+1 , and there are two factors of —1, the first because there is 
a fermion loop and the second because we are calculating the ground-state energy density. 
In summary, our results are 



Fig. la: 5g 



2 {2a + l)!A a+1 



2 a+1 a\ 



Kg.lt: V^+l'T'. (3.6) 

Observe that the two amplitudes above are equal in magnitude but opposite in sign. 
Thus, the final result for the ground-state energy to order 5 appears to be 0. This explicit 
cancellation holds for all values of a; although we expected to differentiate with respect to 
a and set a = 0, there seems to be no need to do so in this case and we obtain directly 

-Aground state 

= + O{5 2 ). (3.7) 

Of course, the cancellation that leads to this result comes from subtracting two logarithmi- 
cally divergent integrals. One may ask whether this cancellation persists if these integrals 
are properly regulated. We examine this question in Sec. [V]. 



B. Calculation of (cfi) 

In contrast, the one-point Green's function arises only from vertices in Eq. (|3.3| ) having 
an odd number of boson lines. The two contributing graphs are shown in Fig. |2|. The 
amplitudes are 

P . ■ x (2a + 2 7 + l)!A a+ T 
Flg ' 2a: ~ m5a 2^(a + 7 )! ' 

Fl g.2b: (3-8) 
& 2 a +^(a + 7-l)! V ; 

Adding these amplitudes, differentiating with respect to a, and setting a = and 7 = |, 
we obtain a nonzero result for the vacuum expectation value of the scalar field: 

(0) =Gx= -i6yj ir A/2 + 0(<5 2 ). (3.9) 
We conclude that this supersymmetric theory has a broken parity symmetry. 
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On the basis of calculations performed in Ref. [11], we believe that (0) remains nonzero 
even when S is a positive integer. One might think that the Lagrangian (1^) is surely 
parity symmetric when 8 is an even integer and one might worry the theory does not exist 
when 5 is an odd integer [because the term — 25+2 appears to be unbounded below]. 
However, neither of these concerns is realized. The reason is that as 5 increases from 
(free field theory) the entire theory must be analytically continued as a function of 5. The 
boundary conditions on the functional integral Z = J Dip exp(— / dxC) representing the 
vacuum persistence function rotate into the complex-0 plane and yield a broken-parity theory 
that exists for all 5. This analytic continuation of boundary conditions is discussed in great 
detail in Ref. 



TTJ. These arguments are based on analysis given in Ref. W% 



C. Calculation of the fermion-boson mass ratio 

Next, we calculate the mass renormalization of the fermion and the boson. These mass 
shifts are obtained by evaluating the negative amputated one-particle-irreducible graphs 
representing to the two-point Green's functions. One graph contributes to the mass shift of 
the fermion to order 6 and two graphs contribute to the mass shift of the boson to order S. 
These graphs are shown in Fig. [| The symmetry number for each graph is indicated in the 
figure. The negative amplitude for the fermion graph is 

Fig .3a: H {2a+ ^f. (3.10) 
6 y 2 a+1 a\ v ' 

The negative amplitudes for the boson graphs are 

Fig. 3b: Sg 



2 (2a + 2)!A Q 



2 a+1 a\ ' 
P J2a + 1)! A a 

Fig. 3c: -5g 2 K . ' ., . 3.11 
6 y 2 a (a-l)! v ; 

Based on our calculational procedure, the next step would be to differentiate each of these 
amplitudes with respect to a and set a = 0. However, this is in fact not even necessary. For 
arbitrary a the fermion mass is 

M fermion = g + 5g [ 2 J +1 ^ + 0(5 2 ) (3.12) 



and the boson mass squared is 



M b 2 oson = 9 2 + ^ (2a 2 + a ^ !AC + 0(5 2 ). (3.13) 
Thus, to order 5 the fermion-boson mass ratio R is unity: 

R= ^~ = 1 + Q{6 2y (314) 
-'"boson 

We conclude from these calculations that while parity symmetry is clearly broken, the 

supersymmetry remains to order 5. In the next section we pursue these calculations to 
second order in 5. 
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IV. SECOND-ORDER CALCULATIONS FOR THE SUPERSYMMETRIC 

LAGRANGIAN 



In this section we extend the calculations of the previous section to second order in 5. 
We begin by expanding the Lagrangian (|1.4|) to order 5 2 : 

£ = \{d<P) 2 + \g 2 ^ 2 + \i$W> + \gH> + \&9^ 

+5g 2 <p 2 ln(i0) + ^SgWHW) + 5 2 g 2 <p 2 ^(i<p)} 2 
+^ 2 ##n(i0)] 2 + \&gW> Hi<f>) + 0(5 3 ). (4.1) 



Next we substitute the identity (|2.6|) and identify a polynomial Lagrangian that in combi- 
nation with the derivative operator Q2.12| ) gives the Green's functions of the theory defined 
by the Lagrangian (|4. 1|) to order 5 2 : 

£ = \(d<P) 2 + l -g 2 <p 2 + + l -g^ 

+ 1(25 + 26 2 - n 2 5 2 a 2 )g 2 4> 2+2a --(25- 26 2 + n 2 6 2 p 2 )g 2 ( p 2+2/3 

+ — (45 + 85a + 25 2 + A5 2 a - Ti 2 5 2 a 2 )g^ 2a 
16 

- — U5 + 85(3 - 25 2 - A5 2 (3 + it 2 5 2 (3 2 )g^<p w 
16 

+ l -m(5a + 5 2 )g 2 (f> 2a+ ^ +1 - -m(5(3 + S^g 2 ^ 2 ^ 1 

+ liir(25a + 5 2 a + 25 2 a 2 )g^<j) 2a+2 ^- 1 - \iv(25p - 5 2 (3 - 25 2 (3 2 )g^(j) 2 ^- 1 . (4.2) 
8 8 

To find this effective Lagrangian requires considerable algebra; one must solve a system of 
eighteen simultaneous equations similar in form to those in Eq. ( j2.13|) . 

In this effective Lagrangian we treat the parameters a, (3, and 7 as positive integers so 
that we can derive a set of Feynman rules for calculating the Green's functions. After these 
Green's functions have been calculated to second order in 5 we apply the derivative operator 
V in Eq. (pT2|) and set a = 0, (3 = 0, and 7 = \. 

The Feynman rules for the Lagrangian ( |4.2| ) are the generalization of the rules in Eq. (|3.3|) 
to second-order in 5: 

boson line: 



p 2 + g 2 ' 



fermion line: 



g-f 



(2a + 2) boson vertex: v x = -(-25 - 25 2 + Tx 2 5 2 a 2 )g 2 (2a + 2)!, 
(2(3 + 2) boson vertex: v 2 = ^(25 - 25 2 + 7c 2 5 2 (3 2 )g 2 (2(3 + 2)!, 
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(2«) 


boson 


and 2 fermion 


vertex: 


^3 = 




(2/5) 


boson 


and 2 fermion 


vertex: 


v A = 






(2a + 


27 + 1) boson 


vertex: 


v 5 = 






(2/5 + 


27 + 1) boson 


vertex: 


ve = 


(2a + 2 7 


-1) 


boson 


and 2 fermion 


vertex: 


v 7 = 


(2/5 + 27 


- 1) 


boson 


and 2 fermion 


vertex: 


v 8 = 



-[-(4 + 8a)5 - (2 + 4a - vr 2 a 2 )S 2 }g(2a)\, 
8 

i[(4 + 8/5)5 - (2 + 4/3 - tt 2 /3 2 )5 2 M2/5)!, 
_(_ fo _^ 2(2a + 27 + 1)!; 

y (5/5 + 5V(2/5 + 27 + l)! ! 
^[-25 - (1 + 2a)6 2 ]g(2a + 2 7 - 1)!, 
^[25 -(1 + 2/3)5^(2/3 + 27-1)!. (4.3) 



A. Calculation of the ground-state energy to second order in 5 

There are thirty graphs that contribute to the ground-state energy density. These graphs 
are organized into seven classes for which the amplitudes combine in a natural way. First, 
we examine the four single- vertex graphs that are constructed from vertices Vi, t>2, t>3, and 
V4. These graphs are shown in Fig. 4(a). The sum of the amplitudes for these graphs is 

Vl A a+1 v 2 A l3+1 v 3 gA a+1 v A gAP +1 

1 ~ ~2 a + 1 {a + l)\ ~~ 2/ 3 + 1 (/3 + l)! + 2 Q a! + 2^\ ' ^ 

When these amplitudes are combined we find that the contribution of order 5 is identically 
0. Only terms of order S 2 survive. Simplifying the expression in Eq. (|4.4| ) gives 



ft 2 n 2 \ a+1 (2n>V 

Ax = 2a+ 3 a! \- Aa + 2 - "M 4 " + !)1 + (« - fl- ( 4 - 5 ) 

Next, we consider the two- vertex graphs constructed from v\ and V3, which are shown 
in Fig. 4(b). Let 21 be the number of boson lines connecting these two vertices. Then, for 
each value of I the sum of the amplitudes for these graphs is 

v 2 A 2a+2 - 21 JdxA 2l (x) vmgA 2 ^ 2 " 21 JdxA 2l (x) 



[(a + l-/)!] 2 2 2a + 3 ~ 2 '(2/)! (a + l-/)!(a-/)!2 2Q + 1 ^(2/)! 
_ v 2 g 2 A 2a+2 - 21 J dx A 21 (x) v 2 A 2a ~ 21 J dx A 21 (x)Tr S(x)S(x) 
[(a-0!] 2 2 2a+1 " 2i (2/)! + [(a-/)!] 2 2 2Q + 2 - 2/ (2/)! ' 

Note that we must sum over all possible values of /. In all of the above terms except the 
last, I ranges from 1 to 00, but in the second v 2 term I ranges from to 00. To simplify the 
Feynman integral in this term representing the boson- and fermion-exchange graph, we use 
the following identity stated in ||: 

f dxA l (x)TrS(x)S(x) = ^ + 2) / dx[A(x)f +2 - I > 0, (4.7) 

J l> I 1 J L I 1 

where S(x) and A(x) are the fermion and boson propagators in coordinate space: 
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S(x) = (g-i$)A(x), 

A{X) = h Koi9lxl) - (48) 

When we simplify the expression A 2 (l) there is a remarkable cancellation that occurs; all 
integrals over A 2l (x) cancel independently of the value of a and we obtain the following 
simple sum: 

oo i->2Z 

A 2 = £ M l) = - 9 ^A^ l( 2a + D!] 2 2— E (a _ 1)!|(a + t (4.9) 

Third, we consider the two- vertex graphs constructed from v 2 and v±, which are shown 
in Fig. 4(c). The sum of the amplitudes for these graphs is identical to A 2 in Eq. (|4.9| ) with 
the replacement a — > (3: 



A, = Ao 



(4.10) 



Fourth, we examine the two-vertex graphs constructed from one of v\ and v 3 and one 
of v 2 and t>4. These graphs are shown in Fig. 4(d). For each value of I the sum of the 
amplitudes for these graphs is 

Vl v 2 A a+ P +2 - 21 JdxA 2l {x) Vl v 4 gA a+ P +2 - 21 JdxA 2l (x) 



+ 



a + 1 -/)!(/?+ l-Z)!2 Q +£+ 2 - 2i (2Z)! (a + 1 - l)\{(3 - l)\2 a +?+ l - 2l {2l)\ 
v 2 v 3 gA a+l3+2 ~ 21 J dx A 21 (x) v 3 v 4 g 2 A a+l3+2 - 21 J dx A 21 (x) 



[a - l)\{(3 + 1 - l)\2 a +P+ 1 ~ 2l {2l)\ (a - 1)1(0 - i)!2«+0- 2t (2Z)! 



v 3 v 4 A a+l3 - 21 J dx A 2l {x)Tr S{x)S{x) 
+ (a -/)!(/? -l)\2 a +P+ 1 ~ 2l {2l)\ ' ( ' ' 

Again, we must sum over all possible values of I; in all of the above terms except the last, 
I ranges from 1 to oo, but in the V3V4 term representing boson- and fermion-exchange, / 
ranges from to 00. Also, we again use the integral identity in Eq. (|4.7|) . As before, when 
we simplify the expression A±(l) there is a remarkable cancellation of the A 2l (x) integrals 
for all values of a and f3. We obtain the following simple sum: 

00 22Z-o-/3-4 

A, = Y.MQ = fPK"»»Va + 1)1(20 + D! g (2 ,_ 1)l(a+1 _,) 1G8+ !_,), ■ <«2) 

Next, we consider the two- vertex graphs constructed from t> 5 and Vy, which are shown in 
Fig. 4(e). For each value of I the sum of the amplitudes for these graphs is 



A 5 (l) 



+ 



v 2 A 2a+2 ^- 21 JdxA 2l+1 (x) 
[(a + 7 - iy] 2 2 2a + 2 ^ 2l (2l + 1)! 

v 5 v 7 gA 2a+2 ^- 21 J dx A 2l+1 (x) 

(a + 7 -/)!(« + 7 - 1 - l)\2 2a + 2 ^- l - 2l (2l + 1)! 
_ v 2 g 2 A 2a+2 ^ 21 JdxA 2l+1 (x) 
~ [(a + 7 - 1 - /)!] 2 2 2a +^- 1 - 2 '(2/ + 1)! 

^20+27-2^ j dx A^(x)Tr g(x)g(x) 
+ [(a + 7- 1 -Z)!] 2 2 2a + 2 ^(2Z + 1)! ' ^ ' 
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Again, we must sum over all possible values of I; in all of the above terms / ranges from 
to oo. Also, in the boson- and fermion-exchange graph in the last term, we again use the 
integral identity in Eq. (|4.7|) . When we simplify the expression A 5 (l) there is no longer any 
cancellation for arbitrary values of a. Thus, the final result still contains an integral over 



A 



21+1, 



We obtain the following sum: 



A 5 = = K 2 g 2 5 2 a 2 A 2a+2 ~<[(2a + 2 7 - 1)!] 2 2" 



2a-27-3 



E 



Y2l 



( , (20![(« + 7 -/)!] 
+ g g 2 2 2l [(2a + 2 7 ) 2 - 1] [(2a + 2 7 ) 2 + 4/ + 1] / dx A 2l+1 (x) } 

1=0 



(21 + l)![(a + 7 -/)!] 2 A 2 ' 



(4.14) 



Sixth, we consider the two-vertex graphs constructed from t> 6 and f§, which are shown in 
Fig. 4(f). The sum of the amplitudes for these graphs is identical to A 5 in Eq. ( |4. 14| ) with 
the replacement 



Ar = A, 



(4.15) 



The seventh and last class of graphs is constructed from one of v 5 and v 7 and one of Vq 
and Vq. These graphs are shown in Fig. 4(g). It is easy to see that the amplitudes for all 
such graphs are proportional to aj3. Thus, anticipating that at the end of the calculation we 
will set both a = and (3 = 0, we need not calculate the amplitude A7 because it will not 
contribute to the ground-state energy density. (Note that before we set a and (3 to 0, we 
must apply the differential operator T> in Eq. ( |2.12| ). This operator does not have a mixed 
derivative term and thus it cannot eliminate both factors of a and f3.) 

The final part of this calculation consists of applying the operator T> in Eq. ( p,12| ) to 



Ax + A 2 + A 3 + A 4 + A 5 + A 6 + A 7 and setting 
lengthy calculation, we obtain the result 



a 



0, /3 — 0, and 7 



After a rather 



7 

ground state ^ t Aj 

3=1 



-5 2 g 2 A 
16 y 



a=0=O, 7=1/2 



2^^)-27T 2 + E 



E 



?r(/) 



1=0 



(/-§)r(/ + i) t!ir(i 



+ 0(<5 3 ). (4.16) 



The sums in this expression may be evaluated easily and we obtain the result 

-Aground state = + 0(# 3 ). 



(4.17) 



B. Calculation of (4>) to second order in S 

The graphs that contribute to the vacuum expectation value of the scalar field arise 
either from one- vertex graphs constructed from the odd vertices, Vc, through t> 8 , or from two 
vertex graphs with one odd vertex and one even vertex. Some of these graphs are shown in 
Fig. |5]. Again, they fall into natural classes. We first consider the five two-vertex graphs in 
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which an even number of bosons are exchanged between the pairs of a vertices (v\, V3) and 
(t> 5 , Vj). For a given I the sum of the five amplitudes is 



21 



v lV7 gA 2a+w - 21 J dx A 



21 



{a + 1 - l)\(a + 7 - l)\(2l)\2 2a +'r+ 1 - 21 
v 3 v 5 gA 2a+2 ^ +1 - 21 JdxA 21 



(a + 1 - l)\(a + 7 - 1 - /)!(2/)!2 2a +^- 2 ' 



(a-/)!(a + 7-/)!(2/)!2 2Q +^ (a - l)\(a + 7 
^A 2 ^- 1 - 2 ' / da; A 2 'Tr S(x)S(x) 
' {a - Z)!(a + 7 - 1 - /)!(2/)!2 2q +7- 2 ' ' 



1 -/)!(2/)!2 



2a+7-l-22 



(4.18) 



There are five more graphs in which an odd number of bosons are exchanged. The sum 
of the amplitudes of these graphs is similar to the result in Eq. (|4.18| ) and we do not give 
it here. Furthermore, there are ten corresponding graphs in which a is replaced by /3 and 
these are constructed from the vertices V2, V4, vq, and v$- When these twenty amplitudes 
are combined and summed over / no dramatic cancellation like that in the calculation of the 
ground-state energy density occurs. Thus, we are left with an infinite sum over integrals of 
the coordinate space propagator A(x). 

Next, we consider the contributions of the twenty graphs, analogous to those above, that 
are constructed from one a vertex and one f3 vertex. That is, we construct all possible 
multi-boson exchange graphs from the vertices (t>i, v 3 ) connected to (t> 6 , v 8 ), and (t> 2 , ^4) 
connected to (t>5, v-j). This calculation simplifies dramatically because only one- and two- 
particle exchange graphs survive when we apply the derivative operator in Eq. ( |2.12| ) and 
set a = (3 = and 7 = 1/2. 

Last, we include the four single- vertex graphs constructed from the vertices W5, Vq, V7, 
and v$. 

When we combine all of these calculations we obtain the following result: 



(0) = ^ttA/2 I -5 + 6 2 



+g 2 J dxA(x) 



'A(x) 
~2A 



1 - In 2 

A s 



ln(^]+V(l) 



+ 



1 + 



A(x) 



A 



In 



1 + 



A fx) 



A 



+0(5 3 ) } . (4.19) 



Note that any positive integer power of the propagator A(a;) in Eq. ( f4.8| ) is integrable. How- 
ever, the value of the propagator at the origin A = A(0) is a divergent quantity. Therefore, 
the function A(x)/A vanishes everywhere except at x — 0, where it is unity. Hence, the 
integral involving this ratio in Eq. ( [4.19| ) exists and vanishes. Thus, our final result for the 
one-point function, which measures the parity symmetry breaking in this theory, is 



(0) 



ttA/2[-5 + 5 2 (1 - In 2) + 0(5% 



(4.20) 



The fact that the theory is supersymmetric reduces the degree of divergence of this result. 
At intermediate stages of the calculation, the coefficient of 8 2 is proportional to A 1//2 In A. 
However, when the boson and fermion contributions are combined, all terms containing In A 
cancel exactly. Thus, the higher-order result is no more divergent than the leading-order 
result. 
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C. Calculation of the fermion-boson mass ratio to second order in 5 



We do not discuss this calculation in detail here because a similar one is explained in 



Ref. [10]. The calculation done here is more elaborate because there are twice as many 



vertices but the necessary calculations are routine. Our result is 



R = Mi ~ = l + 0(6 3 ), (4.21) 

-'"boson 



which is consistent with the theory being supersymmetric. 



V. CONCLUSIONS 

In the Schwinger model of two-dimensional electrodynamics with massless fermions there 
is an anomaly. If the one-fermion-loop contribution to the photon propagator is calculated 
naively one obtains a product of two factors; the first factor vanishes in two-dimensional 
space, and the second factor is a divergent integral. If one is not careful, one gets a quantity 
that is formally 0. However, because the integral is divergent, one must evaluate this prod- 
uct by introducing a regulator; dimensional regulation is an acceptable procedure. As the 
regulator is removed one obtains a finite, nonvanishing result for the anomaly, namely, the 
famous number e 2 /n. In general, one looks for an anomaly when there is a naive cancellation 
involving divergent quantities that must be regulated. The question that is raised in this pa- 
per is, do we have an anomaly in the 5 expansion that breaks supersymmetry? Specifically, 
is there an anomaly associated with the cancellation that gives a vanishing ground-state 
energy density in Eq. ( |3.7| )? 



A. Dimensional regularization of the ground-state energy density 

In the derivation of Eq. ( J3.7|) we combine two numbers that are divergent to obtain 
0. There are several ways to regulate the integral representing A. For example, if we use 
dimensional regulation and evaluate A in 2 — e dimensions, then for small positive e 

27r jo p z + g z 



1 _« r 1 
-9 



4tt Jo 

= A^ 9V (2) 1 i 1 ~ 2, 

= \g~ e . r ,o\ ~T-9- e (e-0+). (5.1) 
4 sm(7re/2J Int 

Furthermore, in (i-dimensional space, the representation of the Dirac matrices has rank 2 d l 2 . 
Thus, the trace of a unit matrix in 2 — e dimensions is 

Trl = 2 (2 " e)/2 ~ 2 - e In 2 (e -> 0+). (5.2) 
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Thus, the coefficient of the second graph amplitude in Eq. Q3.6|) should be multiplied by 
1 - ±eln2. 

We now see that the two graphs do not exactly cancel; rather, the difference in the 
numerical coefficients is of order e. Combining the two graphs in (|3.6| ) now gives 



e(ln2)^ : 



,(2a + 1)!A 
2 a+2 a\ 



a+l 



eln2 

In 2 

4tt 3 / 2, 



5g 2 2 a T [a + - A 



6g 2 T (a+- (2A) 



,a+l 



(5.3) 



If we now differentiate with respect to a and set a = 0, we obtain for the ground-state 
energy density 



E, 



ground state 



In 2 

~8tt 
In 2 

~87 



6g 2 



+ ln(2A) 
— ln(7re) 



+ 0(5 2 ) 
+ 0(5 2 ) 



0+) 



(5.4) 



Because this is a positive number, it suggests that supersymmetry may be broken. 

Of course, dimensional regulation violates supersymmetry. Thus, it is not clear whether 
the nonzero result in Eq. (|5.4|) is correct or is merely an artifact of the regularization scheme 
being used. 



B. Proper-time regularization 

If a supersymmetric regulation exists then, of course, there will be no anomaly. Con- 
versely, if we could establish rigorously that there does not exist any supersymmetric reg- 
ulation of the delta expansion, then there really would be a breaking of supersymmetry. 
We do not know for certain whether a supersymmetric regulation of the 5 expansion exists. 
However, we believe that we have found a relatively simple way to regulate the S expansion 
that is consistent with supersymmetry and thus we believe that there is no anomalous struc- 
ture and that the ground-state energy is truly identically zero. This suggests that while it 
is relatively easy to break parity symmetry, supersymmetry is extremely rigid and is very 
difficult to break. 

Our regulation scheme, which we believe respects supersymmetry invariance, is a variant 
of the proper-time method due to J. Schwinger [TJ|]. (It is well known that the proper- 



time method correctly yields the anomaly in the Schwinger model.) Here, we will use this 
method to define, in an invariant way, the divergent integral A. To make contact with the 
dimensionally regulated result above, let us work in d dimensions: 



(2tt)V + g 2 J (2tt 

Here, s is a proper-time regulator to be taken to at the end of the calculation. We now 
interchange the order of integration in (|5.5|) , and express the momentum integral as the 
product of d one-dimensional integrals: 
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Then, the integral representing A is immediately expressed in terms of a single regulated 
integral: 



A = rdse-* 2 ( 1 \ = j£» fcx - d/2e -, (5 . 7) 



20^/ 2<% d / 2 

If we set d = 2 — e as in the previous subsection, the integral converges when so = 0, and 



we obtain the same result as in flO 



A^r(l). (5.8) 

However, if we want to preserve supersymmetry, we must remain in two dimensions, in which 
case the integral depends logarithmically upon so'- 

1 poo fj T 1 

d = 2: A=— / ^e-^--[ 7 + ln(^ )], (5.9) 

47T Jgs X 47T 

where 7 is Euler's constant. 

The obvious advantage of this regulation scheme is that it treats bosons and fermions on 
an equal footing; the renormalization of the boson and fermion masses is identical. With this 
regulation scheme, the ground-state energy density is zero, as expected. Thus, we believe 
that there is no anomaly in the 5 expansion. 
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FIGURES 



FIG. 1. The two graphs contributing to the ground-state energy to first order in 5. These graphs 
are constructed from vertices having an even number of boson lines. The symmetry numbers are 
shown beside each graph. 

FIG. 2. The two graphs contributing to the one-point Green's function to first order in 5. These 
graphs are constructed from vertices having an odd number of boson lines. The symmetry numbers 
are shown beside each graph. 

FIG. 3. One-particle-irreducible graphs contributing to the mass renormalization of the fermion 
and the boson to order 5. There is one graph for the fermion mass shift and two graphs for the 
boson mass shift. All graphs are constructed from vertices having an even number of boson lines. 
The symmetry numbers for the graphs are shown beside each graph. 

FIG. 4. The thirty graphs contributing to the ground-state energy to second order in 5. These 
graphs are constructed from the eight vertices in Eq. (4.3). We have organized the graphs into sets 
for which the amplitudes combine in a natural way. These sets are labeled (a)-(g). 

FIG. 5. Five of the graphs contributing to {<j>) to second order in 5. These graphs are con- 
structed from the eight vertices in Eq. (4.3). 
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